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Abstract 

We present a consistent analysis of the signal as well as the irreducible 
background for the search of the SM Higgs boson in the ZZ decay channel at 
the LHC. Soft-gluons effects are resummed up to next-to-leading logarithmic 
accuracy, and the results are compared to those obtained with fixed order 
calculations and the MC@NLO event generator. The soft-gluon effects are 
typically modest but should be taken into account when precise predictions 
are demanded. Our results show that the signal over background ratio can 
be significantly enhanced with a cut on the transverse momentum of the 
ZZ pair. We also introduce a fully transverse angular variable that could give 
information about the CP nature of the Higgs boson. 
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1 Introduction 



The elucidation of the mechanism of electroweak symmetry breaking is one of the main 
goals of the LHC physics program. In the Standard Model (SM) and several popular 
extensions such as SUSY, mass generation is triggered by the Higgs mechanism, which 
predicts the existence of (at least) one scalar state, the Higgs boson. The search for the 
Higgs at collider experiments has now being on-going for two decades. The present direct 
lower limit of the Higgs mass in the SM is 114.4 GeV (at 95% CL) [1], while precision 
measurements point to a rather light Higgs, rrih < 200 GeV [2]. 

At the LHC, the main production mechanism will be gg — > H, and if rrih > 180 GeV, 
the Higgs decay into two Z bosons, h — > ZZ, will provide one of the cleanest signatures 
at hadron colliders, i.e., four leptons. Such a final state will allow a very accurate mass 
reconstruction and the best of all possible discovery modes, a sharp peak over a rather 
flat background. At this stage, accurate predictions from theory will be helpful to design 
the best analysis but are not essential to claim a discovery as data alone will provide all 
the necessary information. However, to answer the key questions on the nature of the 
discovered particle, such as its spin, CP nature and couplings, accurate predictions for 
both signal and backgrounds will be required. 

As far as the Higgs signal is concerned, QCD corrections at the next-to-leading order 
(NLO) have been known for some time [3lllj: their effect increases the LO cross section by 
about 80-100%. In recent years, even next-to-next-to-leading order (NNLO) corrections 
have been computed, first for the total cross section [5], and more recently implemented 
in fully exclusive calculations [6l [7]. Note, however, that all the NNLO results use the 
large-mtop approximation, rritop being the mass of the top quark. 

As far as ZZ production is concerned, NLO corrections have been known for some time 
[HI El [ID]- More recent NLO calculations exist that, using the one-loop helicity amphtudes 
of Ref. [TT], fully take into account spin correlations in the Z boson decay [121 fT3] . 

The fixed-order calculations provide a reliable estimate of signal and background cross 
sections and distributions as long as the scales involved in the process are all of the same 
order. When the total transverse momentum of the ZZ pair is much smaller than its 
invariant mass the validity of the fixed-order expansion may be spoiled since the coefficients 
of the perturbative expansion can be enhanced by powers of the large logarithmic terms, 
In^ Mzz/Pt^ ■ In the case of the Higgs signal, the resummation of such contributions has 
been performed up to next-to- next-to-leading-logarithmic (NNLL) accuracy [T^ [T5| [T6]. 

The purpose of the present paper is twofold. We first consider transverse momen- 
tum resummation for ZZ production at the LHC. The resummation of such logarithmic 
contributions was first considered in Ref. [17]. Here we use the resummation formalism of 
Refs. [IHdH] together with the helicity amplitudes of Ref. [TT] (including finite width effects 
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from the Z bosons, but neglecting single-resonant contributions). Contrary to Ref. [T7] we 
fully include the decay of the Z bosons, keeping track of their polarization in the leptonic 
decay. In the large region we use LO perturbation theory {ZZ+1 parton); in the 
region ^ <C Mzz the large logarithmic contributions are resummed to NLL accuracy. 
The present study parallels the one performed in Ref. [I9] in the case of WW production. 
By using these results, we perform a detailed comparison of signal and background cross 
sections and distributions. 

The paper is organized as follows. In Sect. [2] we analyze the impact of transverse 
momentum resummation for ZZ production. In Sect. [3] we compare signal and background 
cross sections and distributions for the search of a Higgs boson of mass rrih = 200 GeV. In 
Sect. IHwe conclude with a summary of our results. 

2 Transverse- momentum resummation for ZZ produc- 
tion 

In this Section we discuss the effect of transverse-momentum resummation for ZZ pro- 
duction at the LHC, and present a comparison to fixed order NLO results obtained with 
MCFM [ISj and to results obtained with MC@NLO p]. 

We consider the process pp ZZ + X e^e~n~^n~ + X and perform the all-order 
resummation of the logarithmically enhanced contributions at small pf^. The implemen- 
tation is completely analogous to the case of WW pair production discussed in Ref. [19] 
and is based on the formalism of Refs. [T5] . We refer the reader to the above papers 
for the technical details. The large logarithmic contributions at small transverse momenta 
of the ZZ pair are resummed up to NLL accuracy. The result is then matched to the fixed 
order LO calculation valid at large pf^, to achieve NLL+LO accuracy. 

We recall that the formalism of Refs. [TH [13] enforces a unitarity constraint such that 
resummation effects vanish when total cross sections are considered. As a consequence, 
at NLL+LO accuracy the integral of our resummed spectra coincides with the total NLO 
cross section if no cuts are applied. 

To compute the ZZ cross section we use MRST2002 NLO parton densities [20] and 
tts evaluated at two-loop order. Our resummed predictions depend on renormalization, 
factorization and resummation scales. The resummation scale parametrizes the arbitrari- 
ness in the resummation procedure, and is set equal to the invariant mass Mzz of the ZZ 
pair. Variations around this central value can give an idea of the size of yet uncalculated 
higher-order logarithmic contributions. Renormalization and factorization scales are set 
to 2Mz- The latter choice allows us to exploit our unitarity constraint and to exactly 
recover the total NLO cross section when no cuts are applied. At NLO we consistently use 
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Figure 1: Comparison of the transverse momentum spectra of the ZZ pair obtained at 
NLL+LO (solid) with MC@NLO (dashes) and NLO results (dots). No cuts are applied. 



l-t-F = f^R = "^Mz as default choice, whereas in MC@NLO /j^p and hr are set to the default 
choice, the average transverse mass of the Z bosons. 

The predictions of resummation are implemented in a partonic Monte Carlo program 
which generates the full 5-body final state (e+e~/i"'"/i~ + 1 parton). Nonetheless, since the 
resummed cross section we use is inclusive over rapidity, we are not able to apply the usual 
rapidity cuts on the leptons. To the purpose of the present work, we do not expect this 
limitation to be essential. 

We start by considering the inclusive cross sections. Our NLL+LO result is 33.76 fb, and 
agrees with the NLO one (33.99 fb) to about 1%. With MC@NLO we obtain 34.60 fb. As 
expected, the MC@NLO cross section is slightly larger because ZZ production is calculated 
in the narrow width approximation, while in the NLO and NLL+LO calculations, finite 
width effects are included. 

In Fig. [T]we show the corresponding p^^ distribution, computed at NLL+LO (solid), 
with MC@NLO (dashed) and at NLO (dots). As is well known, the NLO result diverges to 
+00 as ^ 0, and this divergence is cancelled by the (negative) weight of the first bin, 
due to the virtual contribution. On the contrary, the NLL+LO and MC@NLO results are 
well behaved as ^ — and are very close to each other, showing a peak around ^ ~ 5 
GeV. 

In order to study the perturbative uncertainties affecting our resummed calculation, we 
have varied the renormalization and factorization scale by a factor 2 around the central 
value. We find that the effect of fiR and fip variations is rather small, of the order of ±1%, 
and comparable with the estimated accuracy of our numerical code. Similar effects are 
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Figure 2: Comparison of the transverse momentum spectra of the ZZ pair obtained at 
NLL+LO for different values of the resummation scale Q. No cuts are applied. 



found at NLO. 

The dependence of our NLL+LO results on the resummation scale Q is instead stronger. 
In Fig. [2] we show the NLL+LO prediction for different choices of the resummation scale 
Q. We see that varying the resummation scale the effect on the ^ spectrum is visible and 
amounts to about ±10% at the peak. For lower (higher) values of Q the effect of resum- 
mation is confined to smaller (larger) values of p^^ . Thanks to our unitarity constraint, 
the total rate is instead insensitive to resummation scale variations, within the numerical 
accuracy of our code. 

As in the case of Higgs [U] and WW [123 production, we find that the choice Q = 2Mzz 
gives (slightly) negative cross sections at very large pf^. In order to define a range of 
variations of Q, we prefer to avoid values that give a bad behaviour at large pf^. For 
this reason, in the following, we will consider resummation scale variations in the range 
Mzz/4 <Q< Mzz. 

We now consider the selection cuts designed for the search of a Higgs boson of mass 
rrih = 200 GeV in the e~^e~fi^fi^ channel P2]- The final-state leptons, ordered according 
to decreasing p-r, should fulfil the following thresholds: 

PTi > 22 GeV PT2 > 20 GeV prs > 15 GeV pt4 > 7 GeV, (1) 

the invariant mass of the e~^e~ and the fi^fi~ pairs should be between 

60 GeV < Me+e- ^+^- < 105 GeV (2) 

and the invariant mass of the ZZ pair should fulfil 

190 GeV < Mzz < 210 GeV. (3) 
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Figure 3: Comparison of the transverse momentum spectra of one of the Z at NLL+LO 
(solid) with NLO (dots) and MC@NLO (dashed) results. The lower part of the plot shows 
the NLO and MC@NLO results normalized to NLL+LO. 

With these cuts the NLL+LO result is 5.42 fb, which is about 2% smaller than the NLO 
one (5.51 fb). The cross section from MC@NLO is about 11% larger (6.01 fb). This is 
mainly due to the fact that MC@NLO calculates the cross section in the narrow width 
approximation, and therefore the cuts on the invariant masses of the e^e~ and fi~^fi~' pairs 
are always fulfilled. As in the inclusive case, the effect of scale variations on the rate is 
very small, of the order of ±1%. 

We point out that single-resonant contributions are neglected in our calculation. We 
have used MadGraph/MadEvent |21j to check that these contributions are indeed small 
and found that at LO the effects are smaller than the permille level. Effects from off-shell 
photons pp — i> Z'-f* — > e'^e^^^fi^ are larger. At NLO they decrease the cross section by 
about 1% with the cuts described above, due to negative interference between the Z boson 
and the photon. The shapes of the distributions are, however, not significantly changed. 
Hence, it is safe to neglect these two contributions in the NLL+LO approximation with the 
cuts described above. For selection cuts used in Higgs searches where its mass is smaller 
than the ZZ threshold, the effects from off-shell photons cannot be neglected and have to 
be included. 

In Fig. [3] we show the pj- distribution of one of the Z bosons, computed at NLL+LO, 
NLO and with MC@NLO. Contrary to the ^ spectrum, this distribution is well behaved 
at NLO but the effect of resummation is still visible on its shape. This is evident from the 
lower part of the plot, showing the NLO and MC@NLO result normalized to NLL+LO. 
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Figure 4: Transverse momentum spectra of the leptons: NLL+LO (solid), MC@NLO 
(dashes) and NLO (dots). 



In Fig. m we show the px distributions of the charged leptons, ordered according to 
decreasing pt- Here the NLO prediction is in good agreement with the NLL+LO one. 
MC@NLO, however, predicts shghtly softer leptons. 

The effect of scale variations is still very small for the above distributions. Only in the 
high-p2- tail of the distribution resummation scale variations give a visible effect, being 
of about ±10% at p| ~ Mz- 

In Fig. [5] we consider the distribution in the variable A0t defined as follows. We 
consider the separation between the e~ and the fi' where their momenta are taken in 
the rest frame of their parent Z boson, by neglecting all the longitudinal components. In 
this way the A(j)T is manifestly longitudinally invariant. As will be illustrated later, see 
Sect. [3], this angle is sensitive to the CP nature of a Higgs boson resonance. Due to the 
fully transverse nature of this angle, it can potentially be reconstructed also if only three 
leptons are detected together with missing Ej^. 

We see that the shapes of the NLO and NLL+LO distributions are qualitatively similar. 
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Figure 5: Same as in Fig. [7] but for the A0t distribution. 



Both decrease with increasing separation angle, the NLL+LO prediction shghtly more at 
small angles before it flattens out, while the NLO prediction has a more constant slope. The 
differences are, however, small. The effect of scale variations on the NLO and NLL+LO 
results is again of the order or smaller than 1%. 

We also plotted the prediction for this angle by MC@NLO, although we remind the 
reader that MC@NLO does not include spin correlations in the Z decay. Despite this 
fact, the shape of this distribution is not too different from those obtained at NLO and 
NLL+LO. 



3 Signal and background 

In this Section we perform a consistent comparison of signal and background cross sections 
for the Higgs search in the gg ^ h —>■ ZZ — ^ e~^e~fi~^fi~ channel at the LHC. We consider 
a Higgs boson with mass rrih = 200 GeV and use the numerical program of Refs. [HI [13 
to compute its transverse momentum spectrum. To be consistent with the backgrounqH, 
we work at NLL+LO accuracy and we generate a set of events containing a Higgs boson 
which is then let decay using the MadGraph package [21]. We use the same cuts as in 
Sect. H 

The signal cross section is 7.74 fb. Comparing with the background we get S/B = 1.43. 
In Fig. E] we plot the pt spectra of the leptons for signal and background. We see that the 
spectrum of the leading lepton tends to be slightly harder for the signal, compared to the 

^In our simplified analysis, we consider only the ZZ irreducible background. We neglect other sources 
of reducible background like tt and Zbb which are known to give a much smaller contribution j22j . 
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Figure 6: Leptons px spectra for Higgs signal (dashes) and ZZ background (solid). 



background, whereas the opposite happens for the lepton with the minimum pt- 

In Fig. [7] we plot the AcpT distribution defined in Sect. [2] for signal and background. 
Since this distribution is expected to be sensitive to the CP nature of the Higgs, we also 
consider the case of a pseudo-scalar Higgs boson. As in the case of the scalar, the events 
are generated starting from the transverse momentum spectrum at NLL+LO and then 
letting the Higgs boson decay using the MadGraph package [21]. The computation of the 
spectrum for the pseudoscalar has been done by using a modified version of the numerical 
program of Refs. [HI [13], using the results of Ref. [23^1 . 

From Fig. [7] we see that the shape of the distribution shows remarkable differences in 
the three cases. As shown in Fig. [5l for the background the distribution is rather fiat. 
On the contrary, for the pseudoscalar, the distribution is peaked at central values of A(f)T, 
whereas for the scalar the distribution has a minimum in this region. We conclude that 
this angular variable has a good discriminating potential to assess the CP nature of the 

■''The spectrum for the pseudoscalar at NLL+LO accuracy can be easily obtained by using the fact that 
the real corrections (in the large-mjop approximation) are the same as for the scalar. As such, the only 
difference from the case of the scalar is in the finite part of the virtual corrections [23] , 
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Pt (GeV) 



Figure 8: Comparison of pt spectra of signal and background at NLL+LO, when standard 
cuts are applied. 



Higgs boson. 

We finally consider the possibility to apply an additional cut on the total transverse 
momentum of the four leptons. This idea is inspired by a comparison of the transverse 
momentum spectra of the Higgs boson and of the ZZ pair in Fig. [HI 

We see from Fig. [8] that the Higgs signal is definitely harder than the ZZ background, 
being peaked at ~ 17 GeV. The ZZ background is instead peaked at p-r ~ 5 GeV. As 
such, a cut on the total transverse momentum of the leptons may increase the statistical 
significance. Starting with the standard set of cuts used in the rest of the paper, we 
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Figure 9: The upper plot shows the NLL+LO efficiency as a function of pxcnt for the Higgs 
signal (dashed) and the ZZ background (solid). The lower left plot shows the signal over 
background ratio and the lower right the statistical significance for an integrated luminosity 
oflOft>-\ 

compute the efficiency of the additional cut by defining 

e(pTcut) = f^PT>PTcut/^- (4) 

In Fig. [9] we plot the efficiency as a function of pTcut for the signal and the background. 
We see that the efficiency of this additional cut decreases more rapidly for the background 
than for the signal. As should be expected, the resummation effect is crucial in this case. 
The fixed order NLO efficiencies, not shown in Fig. [9l become unphysically larger than 
unity for small values of prcut- Due to the fact that the efficiency of the background 
decreases more rapidly compared to the signal, the signal over background ratio increases 
with increasing p^cut , as can also be seen from the lower left plot of Fig. O The lower right 
plot shows the statistical significance for an integrated luminosity of 10 fb~^. We observe 
that the statistical significance is maximum when p^cut ~ 15 GeV. 

The latter point, however, requires a word of caution. The predictions presented in 
the present paper are based on resummed calculations obtained in a purely perturbative 
framework. Intrinsic-pr effects are known (see e.g. Ref. [21] and references therein) to af- 
fect transverse-momentum distributions, particularly at small transverse momenta. These 
effects are not taken into account in our calculation. 

As noted in Ref. these non-perturbative effects have the same qualitative impact 
as the inclusion of higher-order logarithmic contributions, i.e., they tend to make the 
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resummed pt distribution harder. The quantitative results shown in Fig. [9] will certainly 
depend on these effects, although the qualitative picture should not change dramatically. 



4 Summary 

Higgs boson production by gluon-gluon fusion, followed by the decay mode h ZZ —>■ 4 
leptons, provides the best discovery channel at the LHC for Higgs masses above 180 GeV. 
For a precise determination of the properties of the Higgs resonance, such as its mass and 
CP nature, detailed theoretical predictions for the signal and backgrounds are necessary. 

In this work we considered a Higgs boson with mass rrih = 200 GeV and performed the 
resummation of multiple soft-gluon emission for the ZZ background. We then compared 
the results with those for the signal in the case of a (pseudo-) scalar Higgs boson. The 
effects from the resummation of soft gluons are modest for observables like the transverse 
momentum of the final state leptons. However, the transverse momentum spectra of the 
Z bosons and of the ZZ pair are sensitive to these effects. 

An angle A(j)T that is sensitive to the CP nature of the Higgs signal is also introduced. 
This angle is defined in a fully transverse way, such that it is longitudinally boost invariant. 
This angle can potentially be reconstructed also if only three leptons are detected together 
with missing Et- 

We also argued that an additional cut on the transverse momentum of the ZZ pair may 
significantly increase the signal over background ratio and the statistical significance. The 
impact of resummation is of course crucial in this case. The above cut could be helpful to 
claim an early discovery or to obtain an easier determination of the nature of the discovered 
particle. 
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